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Abstract 

The conductance of finite systems plays a central role in the scaling theory of localiza- 
tion (Abrahams et al, 1979). Usually it is defined by the Landauer-type formulas, which 
remain open the following questions: (a) exclusion of the contact resistance in the many- 
channel case; (b) correspondence of the Landauer conductance with internal properties 
of the system; (c) relation with the diffusion coefficient D{uj, q) of an infinite system. 
The answers to these questions are obtained below in the framework of two approaches: 
(1) self-consistent theory of localization by VoUhardt and Wolfle, and (2) quantum me- 
chanical analysis based on the shell model. Both approaches lead to the same definition 
for the conductance of a finite system, closely related to the Thouless definition. In the 
framework of the self-consistent theory, the relations of finite-size scaling are derived and 
the Gell-Mann - Low functions /3(g) for space dimensions d = 1, 2, 3 are calculated. 
In contrast to the previous attempt by VoUhardt and Wolfle (1982), the metallic and 
localized phase are considered from the same standpoint, and the conductance of a finite 
system has no singularity at the critical point. In the 2D case, the expansion of /3{g) in 
1 /g coincides with results of the a model approach on the two- loop level and depends 
on the renormalization scheme in higher loops; the use of dimensional regularization for 
transition to dimension d = 2 e looks incompatible with the physical essence of the 
problem. The obtained results are compared with numerical and physical experiments. 
A situation in higher dimensions and the conditions for observation of the localization 
law a ex. for conductivity are discussed. 
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Figure 1: The Thouless scaling construction. The infinite 
system is composed of finite blocks of size L; if only the level 
closest to the given energy E is retained in each block, the 
effective Anderson model arises, with the overlap integral 
and the scattering of site energies Wz^. 



1. INTRODUCTION 

The scaling theory of localization |1| is based on 
consideration of the so called "Thouless number" 
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equal to the conductance Gl = ctlL'^"^ of the cubic 
block of size L {gl is the conductivity and d is a 
dimension of space) in units of /%, or to the ratio 
of parameters Jl and Wl of the effective Anderson 
model, arising in the Thouless scaling construction 
(Fig. 1). Equivalence of two representations in (1) 
follows from the estimate of overlap integrals Jl ~ 
Ti/td through the diffusion time td = Li^ /Dl, esti- 
mate of Wl as the mean level spacing ~' l/vpL'^ 
and the use of the Einstein relation cr^ = e^vpDL 
between the conductivity gl and the diffusion con- 
stant Dl {vf is the density of states at the Fermi 
level). 

The behavior of gL at large L is of the main in- 
terest: if ^ oo, then a system is in the metallic 
phase, since eigenfunctions of blocks are hybridized 
with practically equal weights; if ql — >■ 0, then a sys- 
tem is an Anderson dielectric (hybridization of the 
block eigenstates is practically absent). The block 
of size nL can be composed from vf" blocks of size 
L, so gnL can be recalculated through a given gL as 
gnL = F{gL,n)^ which for n — > 1 can be written in 
the differential form 
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i.e. in the form of the Gell-Mann - Low equation [2] 
The asymptotic behavior of I3{g) 
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Figure 2: (a) To derivation of the Landauer formula (4); (b) 
The difference between Eqs.4 and 5 is determined by the fact, 
that voltage Ues is measured between two reservoirs, while 
voltage Ul between two ideal leads; (c) The many-channel 
scattering matrix. 



follows from the evident relation Gl — o-qoL'^^^ in 
the metallic phase and the estimate Gl ~ exp{~const- 
L} for a dielectric. For d < 2, the /3 function is al- 
ways negative indicating localization of all states. 
For d > 2, it has a root gc, corresponding to the An- 
derson transition point with the power law behavior 
(J oc r" of the conductivity against the distance r to 
the critical point. 

The qualitative considerations of the paper [1] 
stimulated attempts to formulate them in a more 
quantitative form. Conductance of finite systems 
became a subject of a vivid discussion [3]-[Il] (see 
a review article [TS]) resulted in establishing of the 
Landauer approach [5j [12] as an adequate way of de- 
scription. This approach reduces the kinetic problem 
of conductance to the quantum-mechanical scatter- 
ing problem. 

The original Landauer formula for the strictly 
ID (one-channel) conductor follows from the simple 
considerations. If the unit flux of electrons is inci- 
dent from the left to the sample under consideration 
(Fig. 2, a), then it is transmitted with probability T 
and reflected with probability i?=l — r(Tisa 
transmission coefficient). The current through the 
system is proportional to T, while the difference of 
chemical potentials is determined by the difference 
of electron density on the left {1 + R) and on the 
right (T), i.e. l + R-T = 2{1-T). Consequently 



2 



the conductance is proportional to T7(l — T), and 
estimation of the coefficient gives [S]|l| 



Gr = 



T 



2TThl-T 



(4) 



A somewhat different result was obtained by Econo- 
mou and Soukoulis [B| from linear response theory 



Ges = TT^T. 



(5) 



Subsequent investigations established [121 US] ; that 
Eq.4 corresponds to the four-probe, while Eq.5 to 
the two- probe measurement geometry (Fig. 2,b): the 
difference between them is determined by the con- 
tact resistance 27r?i/e^ between the reservoir and the 
ideal conductor 
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Transition from one-dimensional to d-dimensional sa- 
mple requires consideration of the many-channel scat- 
tering matrix shown at Fig. 2,c: the plane wave of 
the unit amplitude incident to the channel i gener- 
ates the transmitted and reflected waves with ampli- 
tudes tij and Tij in the channel j. The multi-channel 
generalization of Eq.5 has a form [31 [71 [TTl [TS] 
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Subtraction of the contact resistance in analogy with 
(6) gives the numerical results [TB] equivalent to the 
Thouless definition [T7], relating the conductance 
with a reaction to boundary conditions^. However, 
multi-channel generalizations of Eq.4 appear to be 
ambiguous [3l 13 [TTl [IS] , and the problem of correct 
exclusion of the contact resistance remains open. 

The previous discussion escapes of the fact that 
the conductance of a finite system is a poorly defined 
quantity. In a strict quantum mechanical descrip- 
tion, a finite system has a discrete spectrum and its 
ground state corresponds to occupation of the lower 
levels (Fig. 3,a)El. If the ground state does not carry 



^ Index L at Gl in Eqs.4 and 6 denotes "the Landauer 
conductance" and not dependence on L. 

^ The Thouless definition provides the equivalence of two 
representations (1). The overlap integral can be estimated 
as the width of the band occurring from the given level in 
the result of the periodic repetition of the block: it is deter- 
mined by the change from periodic to antiperiodic boundary 
conditions. The Thouless definition is physically satisfactory 
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Figure 3; In strictly quantum mechanical description, a fi- 
nite system has no conductance (a) and no resistance (b). 



a current, then a finite conductance is related with 
transitions to excited states, which are separated by 
a finite gap A; such transitions are absent in the 
limit of zero frequency w, and 



ReGL(t^) = 0, 



0. 



(8) 



It is curious, that in the Aharonov-Bohm geome- 
try (Fig. 3,b) the ground state carries a current, if 
the magnetic flux </> through the ring-shaped sample 
is not equal to the integer or semi-integer number 
of quanta ipQ = Tic/e [TJ; in this case, a persistent 
current flows through the system without external 
voltage and the resistance Rl is also zero. 



Rei?L(w) 0, 







(9) 



(contradiction with Eq.8 is avoided due to the pres- 
ence of hnGL{i^))- 

In fact, this problem is well-known: the formulas 
of linear response theory are complemented by a pre- 
scription that the entering them 5 functions should 
be smeared out by the quantity 7, next the thermo- 
dynamic limit i — )■ 00 is taken and only then 7 is set 
to zero. In fact, such procedure transforms the dis- 
crete spectrum into the continuous density of states 
(Fig. 4). For a finite system such a procedure be- 
comes impossible and attenuation 7 should remain 
finite. The question arises on the origin of this at- 
tenuation and its dependence on parameters. 

Attempts of discussing this question were made 
in the papers [IH]~[1I]; based on the "shell model" 
developed in nuclear physics for description of cou- 
pling between a discrete spectrum of the "target" 



but requires consideration of distributions |17| . being hardly 
formulated in terms of average quantities. By this reason it 
practically is not used in analytical theory. 

For simplicity, we have in mind non-interacting electrons 
in the random potential and consider one spin projection. 
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Figure 4: The discrete levels (a) acquire width 7 due to a 
finite lifetime (b). If the limit L — > 00 is taken before the 
limit 7—^0, the extended levels overlap strongly and form 
the continuous density of states (c) . 



and a continuous spectrum of scattered particles [22] . 
Unfortunately, the physics of the problem remain un- 
clear in these papers, because the a model formalism 
was introduced at the early stage of consideration. 
In addition, the bare diffusion coefficient was consid- 
ered as a given constant, while it actually depends 
on the degree of the openness of a system. 

In derivation of the Landauer formulas [51 [71 |H1 
115] the indicated problem is avoided in a following 
manner. The system under consideration is con- 
nected with the ideal leads, which can be taken suffi- 
ciently massive; the spectrum becomes quasi-continu- 
ous and attenuation 7 can be tended to zero. There- 
fore, the Landauer conductance corresponds to the 
composite system " sample-|-external leads" and not 
to the system under consideration. This point is es- 
pecially clear from the fact that the matrix elements 
entering the Kubo formula are determined by inte- 
gration over the region of ideal leads. 

The question arises, in what extent the formulas 
(4-7) reflect the internal properties of the system. 
To illustrate it more clearly, let introduce the poten- 
tial barrier between the sample and ideal leads. If 
the height of the barrier tends to inflnity, then the 
Landauer resistance grows unboundedly, while noth- 
ing occurs with the system itself. Contrary, if the 
height of the barrier tends to zero, then the bound- 
ary resistance disappears but the system is danger- 
ously affected by its environment. 

There is one more question. The inflnite sys- 
tem is fully characterized by the diffusion coefficient 
D{uj,q), which generally possesses the temporal and 
spatial dispersion. The conductance of a finite sys- 



tem is evidently related to D{uj,q) but this relation 
is not clear in the Landauer approach. 

Therefore, the following points remain unclear at 
the present time: 

(a) exclusion of the reservoir contact resistance 
in the many-channel case; 

(b) relation of Eqs.4,5,7 with internal properties 
of the system; 

(c) relation of the Landauer conductance with 
the diffusion coefficient D{uj,q) of an infinite system. 

The answers to these questions are obtained be- 
low in the framework of two approaches: (1) self- 
consistent theory of localization by Vollhardt and 
Wolfle [121 [13' ^^'^ (2) quantum- mechanical analysis 
based on the shell model [T51 [H] . Both approaches 
lead to the same definition for the conductance of 
a finite system, closely related to the Thouless defi- 
nition: it gives a new strong argument in favour of 
the self-consistent theory. Further, in the framework 
of this theory, we calculate the Gell-Mann - Low 
functions [3{g) for the space dimensions d = 1, 2, 3 
(Fig. 5). In contrast to the analogous calculation by 
Vollhardt and Wolfle [H] , the /3 function has no sin- 
gularity in the fixed point gc- The latter is related 
with the fact that the metallic and localized phase 
are considered from the same standpoint, so the con- 
ductance of a finite system has no singularity at the 
critical point: it is in agreement with the general 
principles of the modern theory of critical phenom- 
ena [^ng. 

The present theory has the following structure. 
A finite system is topologically quasi-zero-dimensional 
and its effective dimensionality is less than two. All 
states of this system are formally localized and one 
can introduce the finite correlation length S,qd] it 
satisfies the scaling relation 
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analogous to that for quasi-one-dimensional systems 
[27l|28l|29l[30]; ^ is the correlation length of the infi- 
nite d-dimensional system. The diffusion coefficient 
has a behavior typical for the dielectric phase. 



(11) 



and turns to zero in accordance with (8,9). The 
above statements are valid only for closed systems. 
In open systems the finite diffusion coefficient 
arises, and the following result can be derived for 
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Figure 5: Gell-Mann - Low functions I3{g) for d = I, 
obtained in the present paper. 



2, 3 



the dimensionless conductance 

Con 
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(12) 



Replacement of InL by ln(L/C) in the Gell-Mann - 
Low equation (2) allows to represent gL as a function 
of L/£^. The latter can be determined from (10, 12) 
and allows to reconstruct 13(g)- 

The paper is organized as follows. In Sec. 2 we 
discuss the early attempt of scaling by VoUhardt and 
Wdlfle 24^ and formulate the main difference be- 
tween it and the present paper. In Sec. 3 the corre- 
lation length Cod is discussed and the scaling rela- 
tion (10) is derived for c? < 4. In Sec. 4 we consider 
open systems and derive the result of type (12). The 
same result is derived in Sec. 5 from the shell model; 
its physical sense is clarified and dependence on the 
measurement geometry is discussed. In Sec. 6 the 
length dependence of g^ is presented in the clear 
form and f3 functions for c? = 1, 2, 3 are calculated. 
Their expansion in powers of l/g is compared with 
the results of the a model approach [531 135' ; the use 
of dimensional regularization in a models is found to 
be in conflict with the physical essence of the prob- 
lem. In Sec. 7 the obtained results are compared with 
numerical [3T1[32] and physical ^351 [S3] experiments. 
A situation in higher dimensions d > 4 is discussed in 
Sec. 8. In Sec. 9 we summarize and discuss the con- 
sequences of the present study for the conductance 
distribution, spatial dispersion of the diffusion coef- 
ficient and observation of the localization behavior 
cr{uj) oc —iuj. 



2. SELF-CONSISTENT THEORY AND 
SCALING 

Self-consistent theory of localization by VoUhardt 
and Wolfle theory is based on existence of the diffu- 
sion pole in the irreducible four- leg vertex J7kk'(q), 
entering the Bethe-Salpeter equation and playing 
the role of the scattering probability Wkk' in the 
quantum kinetic equation. Using the estimate in 
the spirit of r-approximation, D cx {U)~^, where 
(...) is averaging over momenta, one can obtain the 
self-consistency equation [23l [24] which can be pre- 
sented in the form '3D10 
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{2t:Y [-^iuj/D{uj)] + q^ 



(13) 

(the limits of integration are written for the mod- 
ulus of q). Here E is the energy of the bandwidth 
order, W is the amplitude of disorder, A is the ul- 
traviolet cut-off. Drain is a characteristic scale of the 
diffusion constant corresponding to the Mott mini- 
mal conductivity. 

The metallic phase is possible for d > 2, when a 
value of the basic integral 
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(14) 



is finite for m = 0. Accepting D = const > for 
a; ^ and specifying r as a distance to a transition, 
one has 

D^D„ 



(15) 



so the exponent of conductivity is unity. In the di- 
electric phase one makes substitution D = —iuj(^'^, 
and Eq.l3 determines the correlation length C; in 
particular, for d > 2 



2 < d < 4 
d > 4 



(16) 



l/(d- 2 
1/2, 

The attempt of using Eq.l3 for derivation of scal- 
ing equations was made in |24) and contains two in- 
gredients. 

* Equation of type (13) can be obtained by approximate 
solution of the Bethe-Salpeter equation |23| or by detailed 
analysis of spectral properties of the quantum collision oper- 
ator 1371 . The possibility to neglect the spatial dispersion of 
the diffusion coefficient was justified in I37| . 
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1. Modification of the Einstein relation. Ac- 
cording to [231, Einstein relation is modified in 
the localized phase due to non-local effects and ac- 
quires the additional exponential factor 



(Tl ^ e vfDl e 
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(17) 



To obtain this result, one considers the change of the 
electron density p{x), induced by the scalar potential 
fix), 

L/2 

p[x) = J a{x~ x')(p{x')dx' , (18) 

-L/2 

where a{x — x') is polarizability 

a{x — x') = —e'^i'F [5{x — x') — 

-{20-'exp{-\x-x'\/0] . (19) 

For a closed system, the diffusion current jdi / / (x) = 
^DLdp{x)/dx at the boundaries of the system x = 
±L/2 is compensated by the electric current je{x) ~ 
<7lE, which allows to determine a^. Producing such 
calculations for 



one has 



p{x) 



e vp 



(p{x) = ipo- Ex, 



E[L/2 + i)e-^''^^ sinh(x/0- 



</5oe"-^/^^ cosh(a;/C) 



(20) 



(21) 



and 



E{L/2 + T (^0 



E{L/2 + ± _L 

2e 
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Accepting ipQ — ztE {L/2 + one obtains 
jei±L/2) = e^i^pDL ( 1 + L/2^ ) e"^/* • E 



(22) 



(23) 



in accordance with (16). It is easy to see that this 
result is related with the unphysical response to the 
constant potential ipo , which is a consequence of ac- 
cepted approximations; absence of self-consistency is 
especially clear for ipo ^ 0, when estimation of ctl 
using je{L/2) and je(— L/2) gives different results. [3 



^ In fact, the kernel a{x, x') should be constructed as a 
binary expansion in eigenfunctions of the diffusion operator; 
for open systems, integration in (18) should be taken over the 
whole space. 



Of course, the correct consideration recovers validity 
of the Einstein relation. In the framework of the self- 
consistent theory, a{u},q) and D{uj,q) are effectively 
independent of q [37], providing the local response 
in the coordinate space; the relation between them 
has a local character and cannot be modified due to 
restriction of the system size. The absence of the 
factor exp{— L/^} is catastrophic for the paper [24) . 
since it fails to obtain the result ql ~ exjp{— const- L} 
in the localized phase. 

2. Modification of the self-consistency equation. 
For a finite system, equation (13) is modified by in- 
troducing the lower cut-off. 



E^ 



Dr 



2-d 



-1/L 



d'^q 1 
{2TrY m? + q^' 



(24) 



and rearranged by subtraction of the same equation 
with L = oo: 



D, 



L'o 



1 2-d 



d'^q 
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(25) 



Since Dao t in the metal and Dqo = in the 
dielectric phase, the diffusion coefficient of a finite 
system acquires a singularity at the critical point. It 
is in conflict with general principles of the modern 
theory of critical phenomena [25, 26 , which allow 
a phase transition only in the thermodynamic limit 
L — >■ oo. 

The present theory is also based on Eq.l3, while 
indicated defects are removed in the following man- 
ner. A finite system is topologically zero-dimensional 
and all its states are formally localized, though the 
effective correlation length foD coincides with ^ only 
in the deep of the localized phase (in the metallic 
regime, ^on > L). As a result, Dqo turns to zero in 
both phases and Eq.25 becomes almost satisfactory. 
The correct result is obtained below, 



E>L — Dm.inK 



-y 
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x| ~L 



(26) 



and differs from (25) by replacement of the integral 
by the discrete sum and concretization of the way 
of cut-off: the latter is provided by the oscillating 
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factor e^'^'^, which effectively restricts summation by 
values |q| < 1/L. These modifications are crucial 
for derivation of the result Dl ~ exp{— L/^} in the 
localized phase. 

The application of the quasi-zero-dimensional con- 
cept is not a pure theoretical construction, but al- 
lows to distinguish the real behavior of open and 
closed systems. Absence of singularities in small sys- 
tems can be also verified experimentally. Therefore, 
the present theory differs from [24] on the level of 
observable consequences. 



3. CORRELATION LENGTH OF 
QUASI-ZERO-DIMENSIONAL SYSTEM 

3.1. Dimensions 2 < d < 4 

A finite system is considered as quasi-zero-dimen- 
sional and its correlation length ^qd can be studied 
in analogy with the quasi-one-dimensional case . 
In a finite system, the basic integral (14) is replaced 
by the discrete sum 



/(to) 
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|q|<A 
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m-'^^OD, (27) 



where allowed values of q has a form 27rs/L, and 
s — (si, S2, ■ • ■ , Sd) is the d-dimensional vector with 
integer components Si = 0, ±1, ±2, . . . We accept 
the periodic boundary conditions in all directions, 
which correspond to the closed system (Sec. 4). The 
term with q = provides the divergency of /(to) 
at TO — 7> and the system is always in the local- 
ized regime. It is convenient to make the following 
decomposition 



/(to) = 



1 1 
1? ^ 



1 

1? 



E 

|q|<A 



1 



1 



1 
1? 
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|q|<A 



= /i(to)+/2(to)+/3(0), (28) 



where we separated the term with q = 0, and the 
rest of the sum is rearranged by addition and sub- 
traction of the analogous sum with to = 0. The 
limit A — > cx) can be taken in the second term I2 (to) 
transforming it to the form L^~'^Ho{mL) (neglecting 
contributions ~ to^A''^^). The third term /3(0) can 



be calculated at i — > 00 by the change of summation 
by integration, while for finite L it has a structure^ 



/3(0) = A'*-^ &0 + 61 



a\d-2 



(29) 



where we accepted a = A ^. Substitution of (28, 29) 
into the self-consistency equation (13) gives 



d-2 



0(to a 



bi+HoimL) 



O 



1 



(30) 



where definition r = /W^ — bo coincides with (15), 
since bo corresponds to the value /(O), calculated 
in the integral approximation. According to (30), 
^OD is a regular function of r. Expressing r through 
the correlation length of the d-dimensional system 
(^-1/"^ ^ |t| = ±t) and omitting terms vanishing at 
a — > 0, one has 



d-2 



= H 



(31) 



{z/2n) 



+bi+ 



(32) 



which is the desired scaling relation (10), consisting 
of two branches {cd are positive coefficients intro- 
duced in [3U]). The asymptotical behavior of H{z) 



l/z\ 

H{z) = { -A{z~ z*) 



z < 1 

z -> z* 

2; > 1 



(33) 



is obtained noticing that H{z) at small z is deter- 
mined by the last term in (32), while for large z 
the sum over s is approximated by the integral; the 
regular expansion is possible near the root z* , corre- 
sponding to the critical point. At arbitrary z the 
sum over s can be calculated numerically, giving 
H{z) for d = 3 shown in Fig. 6. Introducing the 
variables 



one has for dependence y{x) (Fig. 7) 



(34) 



7 



■ 1 H(z) 










d=2 




d=l 




i ' ' 


/ ^ 


/ ^ 




d=3^ 





Figure 6: Functions H{z) for d = 1, 2, 3. 
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Fi gure 7) Behavior of $^q]j/L versus ^/L for d = 1, 2, 3. 
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(35) 



The constant bi in (29) is not universal and depends 
on the way of cut-off: a value bi = —0.226 for a 
spherical cut-off {\q\ < A) is used below for d = 3, 
though in general it should be considered as an ad- 
justable parameter. Corrections to scaling can be 
obtained from (30) and have the same form as for 
quasi-one-dimensional systems: it confirms univer- 
sality of their structure argued in [30] ■ 

3.2. Two-dimensional case 

In two dimensions we have 



^3(0) 



27r a 



bi 



27r a 



and the scaling relation has the form 

L 



27r \L 



H 



(36) 



(37) 



with the previous definition of H{z) (Fig. 6). Using 
the asymptotic results 



H{z) = 



-(l/27r) In^, 



z < 1 
z > 1 



(38) 



1/2 

we have y — x for .t ^ 1 and y ^ [In x] for x ^ 1 
(Fig. 7). Below we use the value bi = 0.1780 ob- 
tained for the spherical cut-off. 

3.3. Dimensions d < 2 

For d < 2, subtraction of the term with m = 
is not necessary and the limit A — )■ cx) can be taken 
immediately in (27). The scaling relation has the 
form 



Hiz) 



1 



i2n) 
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{z/2nr + \s\^ 



(39) 



and consists of one branch, since the function H{z) 
is positive (Fig. 6). Its asymptotic behavior 



H{z) = 



l/z\ 



z < 1 
z > 1 



(40) 



It can be obtained using the o-representation (see Ap- 
pendix) with the cut-off \qi\ < A. 

'^For the cubical cut-off {\qi \ < A) one has 6i = -0.0314. 
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is obtained analogously and and gives y = x ioi x 'ti 
1 and 2/ ~ a;^^-'')/^ for 2: » 1 (Fig. 7). In the ID 
case the function H{z) can be calculated exactly (see 
Eq.118 below). 



4. SITUATION IN OPEN SYSTEMS 

4.1. Difference between open and closed 
systems 

Consider the effective diffusion equation 
Of 



dt 



(41) 



describing the electron distribution /(x, i). In an 
infinite system the operator — has eigenfunctions 
es(x) — e*''= '' and eigenvalues — Qg- In a finite 
system es(x) and As become non-trivial and deter- 
mine the evolution of the initial distribution /o(x) 



/(x,i)=^A,e 



*e,(x) 



As ^ {fa, eg) . 



(42) 

The difference between open and closed systems can 
be formulated on a very abstract level: in the first 
case, the minimal eigenvalue Aq is zero and corre- 
sponds to the constant eigenfunction 

Ao = , eo(x) = const , (43) 
while in the second case 

Aq > , eo(x) 7^ const . (44) 
In the first case one has from (42) at t — > 00 

f{^,t}\,^^= const ^ if 0), (45) 

i.e. distribution f{x,t) tends to the constant limit, 
equal to the average value of /o (x) in the coordinate 
space; consequently, environment does not affect the 
natural process of relaxation, and the total number 
of particles is conserved. In the second case 



f{^,t)\t- 



Aq eo(x) e 



-DXnt 



(46) 



i.e. distribution over x is stabilized, but its am- 
plitude is decreased due to escape of the particles 
through boundaries of the system; by the same rea- 
son their density near boundaries is less than in the 
center. 



Consider examples. Let the ID system is ar- 
ranged in the interval < x < L. For boundary 
conditions of the Bloch type 



f{L) = /(0)e^^ 



(47) 



the allowed values of q has a form Qs = {2t:s + ip)/L 
with integer s, so the system is closed at — 0, and 
maximally open at ip — n, i.e. for periodic and an- 
tiperiodic conditions correspondingly. For the more 
realistic boundary conditions 



/;(0) = «;/(0), f'AL)^-K}{L) 



(48) 



the system is closed for k, = {) due to the absence 
of flow through boundaries; the maximum openness 
is realized at k = 00, i.e. for the zero boundary 
conditions. 

4.2. Attenuation of electron states and 
flniteness of the diffusion coefficient 

In the open system, electrons can escape through 
the boundaries and their eigenstates has a finite life- 
time. It can be proved quite generally, that it pro- 
vides the finite diffusion constant in the static limit. 

Consider the density correlator 

/C.(r-r') = (Gi+Jr,r')G^(r',r)), (49) 

where and are the retarded and advanced 
Green functions. Eq.49 can be rewritten identically 



/C^(r-r') = j de f 

J —00 J — 



duj' ■ 



1 



1 



E + uj - € + iO 

where pg^g+t^{r — r') is the Berezinsky - Gor'kov spec- 
tral density, whose Fourier transform is related with 
polarizability a(a;, (7) [37] 



Pe,e+uiq) = - 



Imae(a;, q) 



(51) 



Using the definitions of kinetic coefficients and an- 
alytic properties of the response functions ^38j , it is 
easy to show that [37 



2'KVp 



(52) 



where D{Lj^q) is the observable diffusion coefficient. 
It should be stressed, that this result has a general 
character and is not restricted by the metallic phase. 
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The finite lifetime of the electron states leads to 
a replacement of infinitesimal quantities ±iO, enter- 
ing definitions of the retarded and advanced Green 
functions, by ±17; analogous changes occur in (50). 
Reproducing the indicated calculations, we come to 
conclusion that the replacement 



—iuj + 27 



(53) 



should be made in (52), both in the term —iuj and 
in D{bj,q). In the localized phase the following com- 
bination remains invariant 



—iu + 27 
(-io; + 27)^2 



(54) 



which has a simple physical sense: attenuation of 
eigenstates was introduced for the permanent eigen- 
functions, so the correlation length ^ characterizing 
the latter is also unchanged. In the static limit 
D{u},q) is replaced by Dl = 27^^, i.e. the finite 
diffusion constant arises. 

4.3. Modification of the self-consistency 
equation 

The result D(u:) -> (Sec. 3) is valid for closed 
systems, being directly related with the existence of 
the allowed value q = 0; the self-consistency equa- 
tion (13) has a form 



= A 



2-d 



-1 _ f 



(55) 

In the open system, the diffusion constant be- 
comes finite, but the correlation length remains un- 
changed: 



(56) 



The labels (c) and (o) indicate the closed and open 
system, which have the different sets of allowed val- 
ues of q in the sum. Taking the difference of (55) 
and (56), one obtains 



which can be considered as a definition of the dif- 
fusion coefficient for a finite system. This definition 




contains essential freedom, since the choice of the 
open and closed system is a subject of agreement. 
For the Bloch boundary conditions^ (47) one can 
accept as etalons the systems with ip — and Lp = ir^ 



9L 




E 

q 



l2 -I- g2 



(58) 

where we came from io gL, accepting D,„i„A2 = 
l/Tivp. We shall refer Eq.58 as the "Thouless defi- 
nition", since gL is determined by the change from 
periodic to antiperiodic boundary conditions^. We 
accept such a change along only one of coordinate 
axes, remaining periodic conditions in other direc- 
tions: according to Sec. 5, it corresponds to the nat- 
ural experimental geometry. 

The definition (58) provides the exponential de- 
crease oi gL in the localized phase (see Appendix) 



9L 



(l-d)/2 



(47r)'*/ 



V 2 / 



(59) 

The origin of the exponential dependence can be ex- 
plained in the following manner. It is well-known 
|39) . that integration of quickly oscillating functions 



— ^ 00 



(60) 



involves the analytic properties of f{x). If f{x) has 
a jump of the n-th derivative at the real axis, then 



-n-l. 



; in particular, the case n = corre- 



sponds usually to integration in finite limits. If f{x) 
is regular at the real axis, then the integration con- 
tour is shifted to the upper half-plane and the inte- 
gral is exponentially small. 



fi^ ^ exp(— const • uj) 



(61) 



The analogous situation takes place in approxima- 



* In general, the factor L~'^ before the sum over q can be 
replaced by the more complicated normalization factor (see 
Eq.108 below), which can give the power corrections in 1/L 
(if treated inaccurately), and destroy the exponent in (59). 
Such problems are absent for eigenfunctions in the form of 
plane waves, corresponding to the Bloch conditions (47). The 
realistic boundary conditions (48) are considered in Sec. 5. 7. 
^ ® Strictly speaking, the original Thouless definition deals 
with the boundary conditions for the electron wave function, 
and not for the effective diffusion problem. Probably, it is the 
most close correspondence that can be established in terms of 
the averaged quantities (see Footnote 2). 
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tion of an integral by a discrete sum 

/oo 



(62) 



S — — 00 



which becomes clear after the use of the Poisson 
summation formula 1401: 



S — — 00 



oo 

E 

k— — oo 



f{x) e'^^^^l^ dx . 

(63) 

The term with fc — corresponds to the integral 
(62), while the main correction to it has an order of 
exp(— const / h) . Two sums over q in Eq.58 are equal 
in the continual approximation; their difference is 
determined by the main effect of discreteness, which 
has the order of exp(— const ■ L) due to h 1/L. 

The definition (58) can be written in the form of 
the scaling relation 



Ht{z 



gL = Ht 
1 



L 

Cod 



E 



(-1) 



2ti 



(z/27r)2 



(64) 



where we introduced a vector t = (ii, ^2, ■ • ■ , ^d), 
whose components ti run integer values 0, ±1, ±2, . . . 
for z = 2, . . . , d and semi-integer values 0, ±1/2, ±1, 
±3/2, ... for i — 1. The function Ht{z) is always 
positive (Fig. 8) and has the asymptotic behavior 



Ht{z) 



(l/vr) {z/2Tif'^~^^'^ 



z < 1 

z > 1 
(65) 




Figure 8: Functions Ht{z), corresponding to the d- 
dimensional "Thouless definition", for ci = 1, 2, 3: in the scale 
of the figure all three curves coincide. 



5.1. Connection of infinite and finite chains 

Consider the model consisting of two chains, up- 
per infinite {N — >■ oo) and lower finite (Fig. 9, a). We 
accept for simplicity that the chains are described 
by the usual Anderson model 



(66) 



though it is irrelevant for the most part of discus- 
sion; the upper chain is supposed to be an ideal con- 
ductor (e„ = 0). As a perturbation, we include the 
overlap integral V between the sites no and toq of 
two chains (Fig. 9, a). The Hamiltonian matrix is 
block-diagonal in the zero approximation, while per- 
turbation creates non-diagonal elements V in no-th 
row and mg-th column and vice versa; so the matrix 
elements of the perturbation operator are 



(67) 



5. APPLICATION OF THE SHELL 
MODEL 

The shell model was developed in nuclear physics 
for description of coupling between the discrete spec- 
trum of the "target" with a continuous spectrum 
of scattered particles [22; lida et al [IH] suggested 
to use it for consideration of the combined system 
" sample-|-external leads". Below we illustrate this 
approach using the simple models of solid state physics 
and then come to consideration of the many-channel 
case.M 



The matrix Dyson equation G 
in components has a form 



Go + GqVG written 



Sec. 5 contains derivation of Eq.64 by the other method 
and can be omitted by the reader interesting only in results. 



Gnn' — + Gj^j^gVGmon' + G j^„^gV G ngii' , (68) 

where Gnn' is the Green function of the perturbed 
system and G'^^^, is the initial Green function cor- 
responding to two independent chains. The index 
n runs values 1, 2, . . . , N corresponding to the first 
chain, and then values iV+l,iV + 2, ...,iV-|-M 
corresponding to the second chain; the sites of the 
latter will be also numerated by the index m. Equa- 
tion (68) is easily solved: setting n = uq and n — mo 

mon' : 



one obtains the closed system for Gnon' and G 
and then its solution is substituted into (68). The 
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Figure 9: (a) Weak connection V between two chains, and (b) the equivalent scheme in the subspace 1. (c) Two weak 
connections between chains, and (d) the corresponding equivalent scheme, (e) Broken upper chain, and (f) effective transitions 
between its parts. 



complete expression for is rather lengthy, so we 
give only its projection on subspace 1 of the upper 
chain 



' momo nono 



(69) 



and subspace 2 of the lower chain 



y2Q0 



(70) 

Investigation of (69, 70) reveals the following quali- 
tative moments. 

1. Effective scatterer. If we are interested only 
in movement along the upper chain, then perturba- 
tion (67) is equivalent to insertion of an impurity 
atom at the point no (Fig. 9,b), with the effective 
Hamiltonian 



Vr. 



nn' - WSnnoSn'no , W = V'Gl„,„„ . (71) 



To prove this result, it sufficient to write down the 
Dyson equation for the perturbation (71) and verify 
that its solution coincides with (69). 

2. Attenuation in the finite system. The initial 
Green function of the lower chain has a form 



^ es(TO)e*(m') 



(72) 



where Cg and es{m) are its eigenvalues and eigenvec- 
tors. In the vicinity of a level the sum is deter- 
mined by one term; its substitution to (70) gives 



es(TO)e*(m') 



(73) 



E-es-V^Gl^^^\eMoW 
For the ideal chain does not depend on n, and 
dk 1 



G\ 





nn 



2tt E- e(fc) + iO 



I{E)-mv{E), (74) 



where v{E) is the density of states at the energy E. 
Since (81) is valid for any level eg, and for small V we 
can neglect the mutual influence of different levels, 
then the effective Green function of the lower chain 
can be written as 



Gn 



sr^ es{m)e l{m') 

2^ F.-f. 



(75) 



E-ls^i'^s 

i.e. the discrete levels acquire a finite attenuation 

7, =^l/2^(e,)|e,(mo)p. (76) 

The difference between 1^ and has no qualitative 
effect and can be neglected for small V . 



3. Effective T matrix. The combination entering 



(69), 



V^G" 



1 _ v^G^ GO 



(77) 
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is in fact the T matrix of scattering [JT]; by defi- 
nition, its substitution into the Born expression in- 
stead of the perturbation V gives the exact scatter- 
ing amphtude. Considering it in the vicinity of a 
level Cs, one finds possibility to write it in the form 



(78) 



It differs from the Born result T = V^G"^^^^ (fol- 
lowing from the concept of the effective scatterer) by 
replacement of G*^ by G, i.e. by taking attenuation 
of states into account. 



5.2. Several bonds between chains 

The simplest generalization of the model contains 
several bonds between chains, connecting the pairs 
of sites uq and mg, ni and mi, and so on (Fig. 9,c). 
In this case the perturbation operator is defined as 



(79) 



and the Dyson equation reads 

(80) 

Let n and n' belong to the upper chain; then Gj^„j . = 
and (80) accepts the form 

Gnn' — G°„/ -I- ^ G^^.VGrmn' ■ (81) 
i 

On the other hand, setting n = in (80) 

Grmn' = y^.G^.^YGn^n' (82) 

j 

and substituting in (81), we obtain 

Gnn' = G^^i -f ^ ^ G'„„. • V G,„.„^ • Gujn' ■ (83) 

If only the upper chain is of interest, then the effec- 
tive perturbation Hamiltonian 

Vnn' = W,j5nn^5n'n, , = V^G^,^^^ 

(84) 

can be used, i.e. the scatterers Wu are introduced 
in the points Ui and the additional overlap integrals 
Wij are included between points rii and rij (Fig. 9,d). 

5.3. Broken upper chain 



Let us remove the portion between uq and ni in 
the upper chain (Fig. 9,e). The above expressions 
formally retain, if the initial Green function G^„/ is 
taken for the broken chain. The equivalent scheme 
shows (Fig. 9,f), that transitions between the left and 
right parts of the upper chain are possible only due 
to overlap integrals Wqi and Wio, so the transmis- 
sion coefficient is proportional to |Wiop. 

5.4. Many-channel case 

Now introduce the model which has the immedi- 
ate interest (Fig. 10, a): the lower chain is replaced 
by a finite d-dimensional system, whose points ■* 
and m^*-* are related with sites n^^ and nj*'' of the 
ideal one-dimensional chains. This model describes 
a situation, when the given d-dimensional system is 
weakly connected with the ideal leads. 

Expressions (79-84) are formally applicable, if 
the index runs values Uq ^ and n^'-* , while the index 



rrii runs values uiq^ and mY' . Therefore, we can im- 
mediately introduce the equivalent scheme (Fig. 10, b), 
according to which the i-th lead on left and the 
j-th lead on right are related by overlap integrals 
Woi, W^i, Wll, Wl{. Comparing with the Lan- 
dauer many-channel scattering matrix (Fig. 2,c), we 
see that the amplitude tjs are determined by the 
quantities W^, 



t 



-2ie 



2ikf 



sink I 



J 



(85) 



where kp is the Fermi momentum. To derive this 
relation, one should write the effective Schroedinger 
equation 

E^pn = Jlpn+1 + Jll^n-1 + ^n1pn + 



+ ^10'5„„<e)^„<.)+W^lf'5„„(»,V'„(.) 



(86) 



and find the solution of the scattering problem: if a 
wave of the unit amplitude is incident to the channel 
s, then the amplitude in j-th channel can be written 
in the form (fc is a wave number) 



^" ~ ' t^je*'=("-"i) , n>ni 



(87) 
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Figure 10: (a) The modelling of the system " sample + external leads". Ideal one-dimensional chains are weakly connected 
with the considered d-dimensional system, (b) Effective transitions between chains corresponding to many-channel scattering 
matrix (see Fig.2,c). 



It leads to the system of equations 



-GoVGoVGoV + ...}\^) 



(91) 



Je 



-ik„ 



J' 



j' 



whose iterations in W^^^ give (85). Substituting (85) 
in the Landauer formula (7), one has 



Gl = 7^4sin2/cF^ VIg" « 
2Trh ^\ ™" 



y4 



J2 



(89) 



This result corresponds to the Born approximation: 
to obtain the complete result, one should find the 
many-channel T matrix. 

5.5. T matrix in the many-channel case 

If $(r) = e*'^ '' is a plane wave and ^'(r) is a solu- 
tion of the scattering problem, then they are related 
by the Lippmann-Schwinger equation |41] 



1^-) = 1$) +60^^!^-) , 
which can be iterated as 

I*) ={\ + GqV + GqVGqV^ 



(90) 



Let us set Go = Gi -f G2, where G2 corresponds 
to the system under consideration, and G\ to ideal 
leads; then perturbation Y relates only G\ and G2, 
while combinations GiV^Gi and G2V^G2 turn to zero. 
Accepting that states |^) and 1$) belong to subspace 
1, we have 

\^)={\ + GyVG2V+ 

+G1VG2VG1VG2V +...}!$). (92) 

By definition, l^*) and |$) are related by the S ma- 
trix, I*) = S\^), and (92) gives 



H-GiyG2y 



1 



l-GiyG2T^ 



(93) 



The T matrix is introduced by the relation y |^) = 
T 1$) m] , so 5 = 1 GoT, reducing to S' = H- GiT 
in the subspace 1; so 



T ^VG2V 



V- 



1 



I-G1VG2V 
1 



-.V . 



E-H2-VG1V ' ^^^^ 

where the relation G2 — (-E — i?2)~^ is used. The 
poles of the T matrix are determined by eigenvalues 
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of the operator H2 + VGiV, which can be found 
perturbatively 



Xs=es + {es\VGiV\es) . 



(95) 



Using the specific form (79) of the matrix elements 
oiV, we have Ag = eg — ijs, where 7^ are determined 
by expression 



7s = TTV'^i^F^\es{mi 



(96) 



which is a natural generalization of (76). This result 
can be also obtained by induction, including connec- 
tions one after another and neglecting their influence 
on Gl,. 

Below we are interested in the limit of small V. 
In this case we can take into account only the quali- 
tative effect related with attenuation, neglecting in- 
fluence of perturbation on eigenfunctions and eigen- 
values. In the Born approximation, the T matrix 
has a form VG2V (see (94)) and has the poles at 
Cg. Substitution of e., — i'js for corresponds to 
replacement of G° by G in expression (89) 



Gr 



2nn 



4sin^A;; 



J2 



r±,r', 



G(r,r') 



\x' —x\—L 



(97) 

where G is deflned analogously to (75). Here we 
introduced the longitudinal (x) and transverse {r±) 
components of vector r, and summation occurs over 
the points of connection with leads. Taking the aver- 
age value of the conductance, we can consider (|Gp) 
as a zero-frequency limit of the density correlator 
(49) and use it in the form (52). Consequently, we 
have related the Landauer conductance with the dif- 
fusion coefficient D{oj,q). Since (97) corresponds to 
the open system, the replacement (53) is implied, 
leading to a finiteness of the diffusion coefficient ■ 
Omitting (here and later) irrelevant constant factors 
we have for the dimensionless conductance 



9l 



r±,r', 



K{r,r') 



1 ^ e*'i ('-''') 



(98) 



5.6. Conductance of a finite system: 
definition 



(a) 





(b) 



Figure 11: The "thin" (a) and "bulk" (b) contacts attached 
to the system. 



Let a finite system has a form of the d-dimensional 

cube connected to external leads, composed of Nc 
ideal one-dimensional chains; one should differ the 
"thin" and "bulk" contacts (Fig. 11). In the first case 
all chains are connected to the spot of a size I ^ 
(where ^/L^ is a characteristic scale where K{r,r') 
essentially changes as a function of the transverse co- 
ordinate rj_), in the second case they are uniformly 
distributed along the side of the cube. For the "thin" 
contacts we can set r± =0 and write 



9L = JlNc^ K{x,x') 



|x-x'|=L 



(99) 



In the metallic state |eg(m)p L~'^ and according 
to (96) attenuation of all states has the same order 
of magnitude 



7 ~ V^vpN^L- 



y2 



(100) 



where A is the level spacing. It is convenient to 
introduce the parameter 



T/2 



(101) 



having a sense of the effective transparency of an in- 
terface. According to (99), ql contains a factor 
in the explicit form and dependence Dl ocy oc kb in 
the diffusion constant, so cc kb- The proportion- 
ality coeflScient can be estimated from the condition 
that for fcft ~ 1 attenuation 7 is of the order A and 
according to scaling theory (see (1,2)) the block of 
size L is in the critical regime, i.e. ql ~ 1 and 



QL^kbL"-^ K{x,x')\ 



\x — x'\=L 



(102) 
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Figure 12; Conductance of a finite system versus degree 
of its openness. 



(a) 





4 













(b) 



This result is valid for kb < 1, when perturbation 
theory is applicable. In the region fc^^l, one expects 
the absence of the kb dependence, since 7 > A and 
the extended levels overlap strongly and form the 
practically constant density of states (Fig. 4) . It is 
easy to see (Fig. 12), that the conductance of the 
maximally open system is obtained from (102) at 
kb - 1.0 

However, the factor kb can be eliminated 
from (102) not only setting kb ^ 1, but also taking 
the derivative at fc;, 0: 



open dgL{kb) 

9l - 



dkb 



= L 



d-2 



■x'\=L 



(103) 

We accept (103) as a definition for the conductance 
of a finite system: physically, it corresponds to the 
extremely open system, but is formulated in terms of 
almost closed systems. Due to the latter, such defini- 
tion reflects the internal properties of the given sys- 
tem, not disturbed by its environment. Simultane- 
ously, it provides the elegant solution of the contact 
resistance problem (Sec. 1): in the small kb limit one 
can use the two-probe formulas (5,7) of Economou- 
Soukoulis type (due to tij — ^ 0) and there is no need 
in the original Landauer formula (4) or its ambigu- 
ous multi-channel generalizations [3l [3 [TTl [13] . The 
allowed values of q in the sum (98) for K{r, r') cor- 
respond to the closed system and include the value 
q = 0, so diverges at to — > and the conductance 
of the ideal system {^qd — 00) appears to be infinite. 
It brightens one of the widely discussed questions 

In the localized phase, transition from (99) to 



Physically, the most adequate estimate of ql for the ex- 
tremely open system corresponds to the plateau at fej > 1 in 
Fig. 12; however, it is not reasonable to use very large values 
of fcj,, since the "plateau" can in fact be a slow ki, dependence 
due to influence of environment on the system. Note, that nu- 
merical modelling |21) usually deals with the limit — f 00. 



Figure 13: A situation in the localized phase. The contacts 
attached to the system arouse metallization of a layer of width 
~ ^ around them (a). If the metallized region is replaced by 
the ideal conductor, then the contacts shift into the deep of 
the system and arouse metallization of the next layer of width 
~ « (b). 



(103) requires more complicated argumentation. The 
estimate (100) for 7 retains for L<^; according to it, 
condition kb ^ 1 corresponds to the critical regime 
for the block of size ^. Let kb ^ I; then the block of 
size is in the metallic state, and hence the layer of 
width ^ ^ around contacts is metallized (Fig. 13, a). 
Let us approximate this metallic region by the ideal 
conductor: then contacts shift into the deep of the 
system and arouse metallization of the next layer of 
the width ^ ^ (Fig. 13,b), and so on. It is easy to 
see that a condition 7 A is valid at all length 
scales till size L. This picture of successive metal- 
lization is valid for kb ^ 1, but retains marginally 
at fcb ~ 1: in this case the blocks of size ^ are in 
the critical regime and the law of their composition 
reduces to stationarity of gL, or Dl ^ JL^^'^, as in 
the metallic phase. For kb ^ 1 metallization does 
not occur and the condition 7 ^ A is valid at all 
length scales. One can see that the estimate (102) 
and the kb dependence (Fig. 12) remain the same as 
in the metallic phase. 

5.7. Equivalence with the "Thouless 
definition" 

The allowed values of q in the sum (98) corre- 
spond to the closed system. Assuming the latter to 
be a system with periodical boundary conditions, it 
is natural to accept its size to be L in the trans- 
verse direction and 2L in the longitudinal direction: 
then for |a; — x'| = L the contacts are arranged at 
the opposite sides of the cylinder, in which the sys- 
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tern is effectively coiled. Considering the ID case for 
simplicity, we have 



1 e*9=^ 



(104) 



qe=2ns/2L 



and noticing that e"?^^ = (—1)* we can separate the 
terms with odd and even s, 



-E 



qs='2ns/L 



qj + 



(105) 



For the Bloch boundary conditions (47) , the allowed 
values are qs = {2-KS + ip)/ L and Eq.l05 contains the 
difference of the terms with = and (p = tt, so 
(103) is equivalent to the "Thouless definition" (58). 

For the more realistic boundary conditions (48), 
the eigenfunctions of the operator —d^/dx^ has a 
form As sm.{qsX + V's) with 

o 

ips = aictan{qg / k) , 



' L + 2k/ {ql + k2) ' 

(106) 

and the allowed values of qs are determined by equa- 
tion 

qsL+2ascia3i{qs/ k) = -KS , s = 1,2,3,... (107) 
The expression for K{x, x') has a form 

^2 siQ(gs^ + V's)sin(g^x^ + '0^) 
K{x,x)-2^As -^^^ , 



and for the closed system (k — ^ 0) reduces to 



qi + m-' 



(108) 



(109) 



-■Ksj L 



We have transformed the product of sines into the 

difference of cosines and extended the summation 
to negative s, using evenness in qg. Noting that 
cos qsL = (—1)*, one can separate odd and even s 
and obtain the result coinciding with (105) apart 
from the irrelevant constant factor. The accepted 
limitation by one dimension is not essential: the 
d-dimensional case differs only by summation over 
transverse components of q, which is the same for 
two terms of the difference (105). 



Below we had in mind the case of the "thin" con- 
tacts (Fig. 11, a). For the "bulk" contacts (Fig. ll,b) 
we have instead (103) 

^o.en^^.-2 1 J2 K(r,r')||._.,|=^ . (110) 



If the one-dimensional chains are connected to each 
site on the plane of the cube, then Nc = L'^~^ and 
summation over r±_ , r'^ removes the transverse com- 
ponents of the vector q, so (110) reduces to the result 
for the ID case. We see that the natural definitions 
for the; conductance of a finite system are exhausted 
by d- dimensional " Thouless definitions" . The intrin- 
sic d-dimensional case is realized for the "thin" con- 
tacts (Fig. 11, a). The effective dimensionality is di- 
minished by the unity if one-dimensional chains are 
connected along the line, which goes through the 
whole plane of the cube. For the "bulk" contacts 
(Fig. ll,b) the effective dimensionality is unity. 

We should notice, that the physical considera- 
tions define to the factor of the order of unity. 
Such uncertainty is natural and related with an arbi- 
trary choice of the unit scale. Only ratios of conduc- 
tances arc relevant, while the choice of the absolute 
scale if a subject of convention. 

We see that one of two scaling relations (10, 12) 
(written as (66) in Sec. 4) can be obtain from the 
pure quantum mechanical consideration without use 
of the self-consistent theory of localization. The sec- 
ond scaling relation can be also studied by other 
methods: in this case, the quantity ^od should be 
defined by Eq.ll, where £>(a;,0) is the diffusion co- 
efficient of the closed system. 

6. DISCUSSION OF SCALING 
EQUATIONS 

According to Sees. 3,4, dependence of ql on L/^ 
is represented in the parametric form 



d-2 



±Cd(^-^j ^H{z), 9l = Ht{z) (111) 

for 2 < < 4, and 

^ln(^|) =if(z) , QL^Hriz) (112) 

for d = 2; in (i < 2 dimensions, representation (111) 
holds with the upper sign and consists of one branch. 
Using asymptotic behavior of II{z) (33, 38, 40) and 
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Figure 14: Conductance versus L/^ for d = 1, 2, 3. 

Ht(z) (65), one obtains for the length dependence 
of 5L (Fig. 14) at d > 2: 

r , ffi » 1 

9L={ g,±B(L/iY-\ gL^g, , 

i (lA)(i/2<)('^-3)/2e-V?, 

(113) 

where = Ht{z*), B = CdH^{z*)/ H'{z*). For d < 
2 the resuhs for small and large g^ are formally the 
same, but the critical point gc is absent. For d = 2, 
the result in the metallic phase gL — (I/Stt) ln(^/L) 
can be represented as a logarithmic correction to the 
Drude conductance go, gL = go — (I/Stt) ln(L/a) [1], 
with the asymptotics (113) for g^ ^ 1. 

The Gell-Mann ~ Low function f3{g) is deter- 
mined by the derivative dlng/dlriL (see Eq.2) and 
can be written in the parametric form 



9 = Ht{z), 
for 0? 7^ 2, and 
g = Hriz) , 



/3(5) = (rf-2) 



H{z)H^{z) 
Ht{z)H'{z) 



1 H'j,{z) 
2^ H'{z)Ht{z) 



(114) 



(115) 



for d = 2. Since Ht{z) is positive, while H^{z) 
and H'{z) are negative (Figs. 6, 8), the (3 function is 
negative for d = 2, and has a root for d > 2 due 
to the root of H{z). The calculated /3 functions for 
d = 1, 2, 3 are shown in Fig. 5. 

The expansions of H(z) and Ht{z) in powers of 

has a following form for d < 2 

H{z) — l/z^ + flo — a2z'^ + a4Z^ — a^z^ + . . . , 
Ht{z) = l/z'^ + ao-a2z'^ + a4z'^-dGZ^ + .. . , (116) 



a2r. 



E 

s=iO 



i2n\s\ 



2ri+2 



? (27r|t 



(117) 

and vectors s and t are the same as in (28) and (64). 
In the case d > 2, the coefficient qq is replaced by 
the parameter bi, introduced in (29) and (36). In the 
case d = 1 , the coefficients a2n and d2n are expressed 
in terms of the Riemann ^ function or the Bernoulli 
numbers [JD] and can be obtained from the closed 
expressions 

^^""^ " 2z tanh(z/2) ' ^^^^^ " z sinhz ' 

(118) 

following from (39) and (64) with the use of the Pois- 
son summation formula. Using (116), one can find 
the expansion of (3(g) in powers of 1/g 



P{g) = (d-2) + 



1 



Zirg 



(d-2)(Qo-ao) 
9 

02 — 02 



2TTg3 



d^2 (119) 



(120) 



The latter result can be compared with the expan- 
sion obtained in the a model approach 1331 [31] 



'2f + 0-t^ + Q-t^ - 12C(3)t^ + . . . , (121) 



which is written in terms of the variable t ^ l/g. 
Recalculating (120) to the same form, one has 

/3(t) = -2t^ + 0-t^ + 327r2(a2-d2)t^ + . . . . (122) 

The first two coefficients coincide with (121), while 
the third one depends on details of the g^ defini- 
tion, since the parameter 0,2 is different for the "thin" 
and "bulk" contacts (Fig. 11). Such situation is well- 
known in the quantum field theory [4311441 145] . where 
the structure of expansion for the /? function is the 
same as (121,122): the first two coefficients are in- 
variant, while the rest depend on the renormalization 
scheme. Transformation from one scheme to another 
corresponds to the change of variables i — f{t), re- 
lating two different definitions of the charge; expan- 
sion of f(t) in the series and the proper choice of 
the coefficients allows to transform (122) into (121) 
[m US]. The function f{t) is well-defined in per- 
turbation theory but can be singular at i ~ 1, in- 
dicating that one of two schemes is surely defective 



These results can be used for numerical calculations in 
higher dimensions, in order to produce the analytic summa- 
tion along one of coordinate axes. 



18 



[IS]. Since (122) corresponds to the physical def- 
inition oi gL, such problems can refer only to the 
expansion (121). In the framework of perturbation 
theory Eqs.121 and 122 are completely equivalent. 

Such equivalence is destroyed in the space dimen- 
sion d = 2 + e. The dimensional regularization used 
in the a models corresponds to the /? function of the 
form 

/32+e(5) = e + /52(ff), (123) 

i.e. the d dependence is present only in the first term 
of l/g expansion. The exponent v is determined by 
the derivative of P{g) at the fixed point, and the 
corresponding result [32 

i. = lA-(9/4)C(3)e2 + ... (124) 

is in conflict with (16). This fact is usually consid- 
ered as a proof that the Vollhardt and Wolfle theory 
cannot be exact. 

However, another interpretation is possible. Let 
us assume (in accordance with [HI [57]), that the 
Vollhardt and Wolfle theory is correct: then the 
physical reality consists in existence of the exact re- 
sult v — 1/e for d = 2 + e and the non-trivial (3 
function for d ~ 2. The latter is related with the 
physical essence of the problem: the logarithmic be- 
havior at g ^ 1 (see Eq.3) makes impossible for 1/g 
expansion to be truncated at finite number of terms. 
Such physical reality is incompatible with the for- 
malism of dimensional regularization: according to 
(123), the exact result = 1/e is possible only for 
the trivial function P2{g) — A/ g. Description of re- 
ality with such formalism should lead to unsolvable 
problems. Exactly such situation takes place in the 
modern theory: the Anderson transition problem re- 
duces to the a model in a certain approximation but 
the corresponding renormalization group is unstable 
to high gradient terms [37] . It is interesting to to 
carry out renormalization of a models without the 
use of dimensional regularization: there are indica- 
tions that in this case the high-gradient catastrophe 
is absent (see discussion of HH] in the paper [47)'). 

The latter asymptotics in (113) can be compared 
with the exact results for the ID case |49] 

{g) = ^{aL/n)-^'^ e~"^^'\ 

exp(ln.g) =4e""^, (125) 



where a oc W'^ for weak disorder. The effective cor- 
relation length in the dependence exp(— L/^) is sen- 
sitive to details of the averaging procedure and deter- 
mined only by the order-of-magnitude, in correspon- 
dence with its physical sense. The same uncertainty 
exists in the present theory, where ^ is defined by the 
relation D{uj) = {—ioj)^"^ with the ambiguous choice 
of the absolute scale for D{uj). Within such uncer- 
tainty, there is no sense to discuss the precise form of 
the pre-exponential dependence which corresponds 
to redefinition of ^ by the factor 1 -I- 0(lnL/L). 

7. COMPARISON WITH NUMERICAL 
AND PHYSICAL EXPERIMENTS 

Behavior oi gL versus L/(^ for d = 3 is com- 
pared (Fig. 15, a) with numerical results by Zhareke- 
shev where gL was estimated as "acceleration" 
of levels Ks = d^es/dip^ at = ((^ is the pa- 
rameter in the boundary condition (47)), which is a 
possible variant of the Thouless definition [17] . The 
agreement is satisfactory in the metallic state and 
the vicinity of transition, while there are expectable 
deviations in the localized phase: they are related 
with the fact that theoretical results correspond to 
ln(5), while numerical to (In 5). According to (125), 
the difference of two situations corresponds to redef- 
inition of ^ by the constant factor, which reduces to 
the parallel shift in the logarithmic scale of Fig. 15, a. 

Comparison of the same dependence with the 
physical experiment [53] is possible under assump- 
tion that L is replaced by the length Lin oc T~", 
characterizing the inelastic processes. Unfortunately, 
there are no grounds for a to be the same in the 
metallic and localized phase, and it can have a slow 
drift as a function of disorder. In Fig. 15, b it was 
suggested that a is the piecewise constant quantity, 
taking different values in the metal and insulator. 
Such assumption does not strongly affect the results 
in the critical region where the length dependence of 
(7l is rather slow. The latter region is poorly pre- 
sented in Fig. 15, b, and in fact it illustrates a situa- 
tion not very close to the critical point. On the other 
hand, the critical behavior obtained in |53| , is excel- 
lently described by the Vollhardt and Wolfie theory: 
the values s = 1.0 ± 0.1 for the conductivity expo- 
nent and z = 2.94 ± 0.3 for the dynamical exponent 
agree with the theoretical results s = 1 and z = 3. 

The length dependence oi gL for c? = 2 can be 
compared (Fig. 16, a) with the numerical data by Mar- 
kos [3T]. There is a good agreement in the region 
W > 2, while for weak disorder numerical results 
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Figure 15: (a) Comparison of theoretical scaling curves for d = 3 with numerical results by Zharekeshev |32| . The general 
form of the curves was determined without adjustable parameters, and only parallel shifts along two axes were made (the 
same for two branches), corresponding to the choice of the absolute scales for gj^ and ^. (b) Comparison of the same curves 
with the experimental results for Si-P 1531 under assumption L oc T~°', where a was chosen independently for two branches. 



g 




Figure 16: (a) Comparison of theoretical scaling curves for d = 2 with numerical results by Markos 31. Fig. 37]. The form 
of the curves was determined without adjustable parameters, only parallel shifts along two axes were made, (b) Comparison 
of the theoretical /3 function for d = 2 (Fig. 5) with the empirical /3 function extracted from experiment 35 under assumption 
L oc T~°' . The open and dark symbols corresponds to the freshly cleaved surface of Ge and to bicrystals of Ge. 
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display a strong violation of scaling: it can be re- 
lated with the protracted ballistic regime [3T] or ex- 
istence of the 2D metal-insulator transition. The 
latter, according to [55], occurs in roughly half cases 
and belongs to the Kosterlitz-Thouless type, com- 
patible with scaling theory [I]. From viewpoint of 
the general analysis [37], a situation at o? = 2 is 
special and probably reduces to the VoUhardt and 
Wdlfle theory not in all cases. 

The theoretical [3 function for d = 2 (Fig. 5) can 
be compared with the empirical /3 function (Fig. 16, b) 
obtained by Zavaritskaya [35j under assumption L oc 
T^". The use of the constant a allows to describe 
a situation both for large and small g. Some devia- 
tions are present only in the region g ^ 0.1, where 
the experimental results are also ambiguous. 

8. SITUATION IN HIGHER DIMENSIONS 

8.1. Dimensions d > 4 

For d > 4 we have for l2{m) in (28) 

l2{m) ^ -m^ {cdA'^"^ + 0{m'^-\L'^-'^)} , (126) 

i.e. analytical calculation is possible (in the main 
approximation) for arbitrary values of m and L~^. 
Indeed, for m>L^^ the sum can be estimated by the 
integral, which converges at the lower limit already 
for rn = 0, so finiteness of m gives only small correc- 
tions. In the case m < , the main effect of finite 
L is related with absence of the term with q = 0, 
which can be estimated as restriction jql > in 
the integral approximation. 

Substitution to (13) reveals the possibility to ne- 
glect &i and leads to the scaling relation 



X'' 



y — 2 



(127) 



in variables 



y ■ 



ioD /ax('i-4)/4 



L \l) 



!(!) 



(rf-4)/4 



(128) 

where we redefined the scales of ^oz> and ^, in or- 
der to obtain the unit coefficients in (127). Scaling 
relations (127, 128) contain the atomic scale a due 
to nonrenormalizability of theory [3D]. The critical 
point corresponds to ?/ = 1, so that 



$0D 

L 



(d-4)/4 



(129) 



and the small z asymptotics can be used for Ht{z) 
(see Eq.65). In the vicinity of transition we can re- 
place £,i^dIL by ^Jgi, and equations (127, 128) deter- 
mine the length dependence of gL] in particular, at 
the transition point 



9L 



^\ (d-4)/2 

a) 



0. 



(130) 



The physical sense of this result is clarified by the 
fact, that it can be obtained from the self-consistency 
equation for an infinite system (following from rela- 
tions of Sec. 2) 



D{uj) ^ At + B 



Dicu) 



l/2u 



(131) 



if one replace D{lu) — )• Dl, —iuo — > 7 and accept 
7 (for r = 0) to be of the order of level spacing 
A (X L~'^. For d < 4 it gives g^ — const, while for 
d > 4 it reduces to (130). 

8.2. Four-dimensional case 

For d = 4 we have analogously 

{-C4,m^\ii— + 0{\) , mL>l 
m 
-C4TO^ln(AL) + 0(l), mi<l 

(132) 

and two results differ by In(mL), which in the actual 
region reduces to the double logarithmic quantity. 
Neglecting such quantities, we can obtain the scaling 
relation (127) in variables 



y= — iHL/a)] 



(133) 

In the vicinity of transition we can replace £,od/ L by 
^/gT, and obtain 

gL-ML/a)]^''-'^'\ r = 0. (134) 

As was explained in [30], Eqs.127, 133 allow to pro- 
duce the usual constructions of scaling curves, if the 
quantity y is considered as a function of the "mod- 
ified length" ^(L) = L [ln(L/a)]~^^^; then a change 
of the scale for ^[L) allows to reduce all dependen- 
cies for T > and t < to two universal curves. It 
should be emphasized, that the critical point cannot 
be determined by the condition g^ = const. 
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9. CONCLUSION 



possibility of constructions 



The present paper continues the hne initiated by 
the previous publication [30] : since there are serious 
indications [371 US that the VoUhardt and Wolfle 
theory predicts the correct critical behavior, it is de- 
sirable to derive its consequences (as many as pos- 
sible) and compare them with the numerical and 
physical experiments on the level of raw data. Such 
approach already proved its value [5D]: the results 
V = 1.3 1.6, obtained usually for d = 3 in numeri- 
cal papers, can be explained by the fact that depen- 
dence L + Lq with Lo > is interpreted as L^^'^ with 
v > 1, while the raw data are excellently compati- 
ble with the self-consistent theory. The finite-size 
scaling relations for the conductance and Gell-Mann 
- Low functions 13(g) obtained in the present paper 
are also in a good agreement with numerical and 
physical experiments. 

In the present paper, we have elaborated a new 
definition for the conductance of finite systems, bright- 
ening the questions formulated in Sec.l. It appears, 
that both self-consistent theory of localization and 
the quantum-mechanical analysis based on the shell 
model lead to the same definition, closely related 
with definition by Thouless. It gives one more seri- 
ous argument in favour of the VoUhardt and Wolfle 
theory. Expansion of the /3 function in 1/g shows 
that there are no contradictions on the perturba- 
tive level with the results of the a model approach 
in two dimensions. Further, in the case of valid- 
ity of self-consistent theory, the formalism of dimen- 
sional regularization is incompatible with the physi- 
cal essence of the problem. Probably, it is the reason 
both for the high-gradient catastrophe, and contra- 
diction with the VoUhardt and Wolfle theory in the 
space dimension d — 2 + e. 

The new definition will probably resolve the prob- 
lem of pathological singularities in the conductance 
distribution [31j . which cannot exist in finite sys- 
tems. Their observation in numerical studies [31^ is 
probably explained by the fact the considered sys- 
tem was not sufficiently isolated from environment 
and the thermodynamic limit L — >■ c» was effectively 
taken along one of the coordinate axes. 

The above approach suggests the simple argu- 
mentation on the spatial dispersion of the diffusion 
coefficient. It is easily proved for the localized phase 
[37) . that D{ijj,q) has a regular expansion in q^. 
However, it does not exclude the appearance of non- 
integer powers of q at the critical point [SO] due to 



D{io,q) 



1 



(135) 



becoming singular in the ^ — ?■ oo limit. In a finite 
system the role of ^ is played by ^or>7 which is reg- 
ular at the transition point; with such replacement, 
Eq.l35 is valid in the metallic phase. However, the 
absence of such dispersion in the metallic regime is 
easily established from the kinetic equation. Accord- 
ing to ^37 , the following result is valid instead (135) 

D{u, q) = (-tw)e' (1 + diq^ + dag^ + . . .) (136) 

(with di independent of r) , which reveals no patholo- 
gies for replacement of ^ by ^od- It should be noted 
that Wegner's exact result D{ijj, 0) ^ u}^d.-2)/d 
the critical point [FT, following from (131) for t — 0, 
cannot be obtained in the case of essential spatial 
dispersion of D{uj,q). It is interesting, that the re- 
cent experiments on the spreading of wave packet 
[5^ are in agreement with the self-consistent the- 
ory and give no evidence of the anomalous spatial 
dispersion. 

The localization law for conductivity cr(a;) cx —ioj 
was predicted almost 40 years ago [53] but has never 
been observed experimentally. The above analysis 
clarifies that its observation is possible in the closed 
systems under approximately the same conditions as 
for existence of the persistent current in the Aharonov- 
Bohm geometry (Fig. 3) [54 1 155 1 156 ] . 



Appendix. Asymptotics of for mL — > oo 

Consider the sum (27) for the Bloch boundary 
conditions in all directions 



1 



m'^ +ql + 



and introduce the so called a-representation 
1 



+ <f Jo 



(Al) 

dae-<^'+"'') . iA.2) 



Then 



I{m) = / dae-"™' TT Sj( 
■Jo 

OO 

S,{a) = L-^ E 



a 



(^.3) 
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and the use of the Poisson summation formula 
transforms Sj{a) to 



1 



/Ana 



3^3 IS 



kj — — oo 

Then (A. 3) takes a form 
da 



I{m) 



(47ra)'*/2 



(A.5) 

where a vector fc = (/ci, . . . , kd) is introduced and 
k ■ if = kjifj. The term with fc = is calculated 
exactly and corresponds to the continual approxima- 
tion. The main effect of discreteness is determined 
by the terms with |fc| = 1, which can be calculated 
for mL ^ 1 in the saddle-point approximation. Re- 
maining only these terms, one has 



/(m) 



mL 



(47r)d/2 

(l-d)/2 



-mL 



d 



2 cos Lpj 



(A.6) 



Taking the difference of two such expressions with 

tfi=Q and ipi = tt, 



(l-d)/2 



(4^)d/2"* ^ 2 

we come to Eq.59. 



-mL 



(A.7) 
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